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I .  INTRODUCTION 

The  discovery  of  surface  enhanced  Raman  spectroscopy^  has 

Inspired  a  large  body  of  work  which  has  been  siimmarized  in 

2—8 

several  recent  review  articles.  While  much  has  been  learned 

about  various  enhancement  sources  there  is  still  disagreement 
regarding  the  proportion  in  which  different  factors  contribute  to 
the  observed  Raman  cross  section.  Even  though  it  is  difficult  to 
Include  all  the  opinions  in  a  broad  classification,  it  is  fair  to 
say  that  the  main  debate  is  centered  around  the  relative 
importance  of  "chemical"  versus  "electromagnetic"  (EM) 
enhancement  factors . 

The  two  kinds  of  effects  can  be  defined  semi-quantitatively 
by  examining  various  terms  appearing  in  the  expression  for  the 
intensity  of  Raman  scattering  by  a  molecule  located  near  a  solid 

9 

surface.  A  source  of  light  of  frequency  w  producing  an  electric 
field  E^{r,w)  in  vacuum,  will  generate  a  different  field  Ej^(r;w) 

« r ;«) »  if  a  solid  surface  is  present.  When  a  molecule 

is  placed  at  r  it  will  interact  with  and  acquire  a  Raman 

4  5  ^ 

dipole  '  given  by 

*♦ 

=  Q  •  R(r;«) (I.l) 

Here  Q  is  the  amplitude  of  the  normal  mode  of  interest  and  oc  is 
the  polarizability  tensor  of  the  molecule.  The  electric  field  ^ 
of  the  radiation  emitted  by  this  dipole  can  be  written  as 

f(r  =*  G*li  »  Q^''^*R(r;«)  •1^(7 ;w)  (1.2) 

-►  ♦♦ 

where  r'  is  the  position  of  the  detector  and  G  is  the 

4—5 

electromagnetic  Green's  function  for  a  dipole  oscillating  near 

a  surface.  The  Raman  intensity,  I,  corresponding  to  Eq.  (2)  is 

^  2 

proportional  to  | E( r ' : | ' 

Equation  (2)  permits  a  simple  classification  of  the  total 
enhancement  into  terms  involving  different  mechanisms:  the 
enhancement  of  G  and  R  is  usually  electromagnetic,  while  that  of 


2 


3a/3Q  is  "chemical".  Large  enhancements  of  G  and  R  can  be 

achieved^*^  by  using  surface  shapes  and  materials  which  permit  the 

incident  light  and  the  emitting  dipole  to  excite  the 

electromagnetic  resonances  (e.g.  plasmons)  of  the  solid. 

Moderate  of f -resonance  enhancements  of  G  and  R  can  also  be 

1  Oc 

obtained  by  using  high  surface  curvatures  ,  or  the  mutual 

lOe 

polarization  of  two  solids  located  next  to  each  other.  Many 

examples  of  this  kind  of  enhancement  are  described  in  refs. 

2', 4,6-8  and  especially  in  ref.  5.  An  additional  electromagnetic 

effect  is  caused  by  the  coupling  of  the  emitting  dipole  to  its 

image  field,  which  shortens  the  lifetime  of  the  emitting 

state® and  depresses  the  intensity  of  the  resonant  Raman 
12 

signal . 

The  magnitude  of  the  electromagnetic  effects  is  estimated 

by  solving  the  “phenomenological"  Maxwell  equations  under  a 

number  of  assumptions.  Unfortunately  these  assvimptions  are 
13  14 

inaccurate  '  and  their  use  leads  to  errors  whose  magnitude 
cannot  be  adequately  estimated.  This,  together  with  the  fact 
that  in  working  with  rough  surfaces  it  is  extremely  difficult  to 
control  surface  structure  {which  strongly  influences  the  EM 
effects),  makes  it  difficult  to  assess  liow  much  of  the  observed 
enhancement  is  due  to  electromagnetic  effects  and  how  much  to 

5 

chemical  ones.  Nevertheless  a  broad  survey  of  various 
spectroscopic  probes  on  a  great  variety  of  surfaces  shows  that 
the  electromagnetic  theory  is  in  qualitative  agreement  with  many 
experimental  observations.  Furthermore,  some  of  the  most  careful 
studies  in  the  SERS  literature,  carried  out  by  Murray  and 
Bodoff^®,  Mullins  and  Campion^®  seem  to  indicate  that  the 
chemical  effects  are  small. 

However  not  all  molecules  are  alike,  and  a  large  number  of 

studies  have  pointed  out  behavior  at  variance  with  the 

17 

qualitative  predictions  of  the  EM  theory  The  extent  to  which 

these  studies  are  conclusive  is  still  debated,  and  a  large  r amber 

17  —  28 

of  theoretical  models  have  been  used  to  suggest  possible 

sources  of  chemical  enhancement,  and  to  estimate  its  magnitude. 


The  majority  of  these  proposals  are  summarized  by  Fig.  1 
which  schematically  represents  the  Raman  process  for  an  adsorbed 
diatomic  molecule.  The  incident  photon  of  frequency  u  excites  an 
electron  from  the  ground  state  |0>  into  the  state  |i>;  this 
electron  moves  through  the  solid,  excites  the  "local  phonon" 

(e.g.  the  diatomic's  vibrational  stretching  mode)  and  is 
scattered  into  the  state  |j>;  subsequently  the  electron  decays 
into  the  ground  state  and  emits  a  Raman  photon.  The  total 
emission  is  obtained  by  adding  up  all  such  processes  occurring  at 
all  the  points  A,  B  and  C  in  the  system,  for  all  the  electrons 
that  can  participate  without  violating  the  conservation  laws  or 
the  Pauli  principle. 

To  say  that  an  accurate  calculation  of  this  process  is 
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extremely  difficult  is  an  understatement.  '  Thus,  the 
existing  models,  including  the  one  to  be  presented  here,  have  had 
to  resort  to  severe  approximations.  This  makes  them  conceptually 
useful,  but  limits  their  predictive  power  and  reliability. 

The  diagram  in  Fig.  1  allows  us  to  classify  these  models 

into  two  groups.  One  group  relies  on  the  fact  that  the  electron 

lines  (i.e.,  the  one-electron  Green's  functions)  have  energy 

denominators  which  lead  to  resonant  Raman  signals  when  the 

incident  laser  or  the  emitted  light  match  an  electronic 

transition  in  the  molecule-solid  system.  Since  the  excited 

electron  must  interact  with  the  local  phonon  a  substantial 

enhancement  of  the  Raman  signal  can  come  only  from  transitions 

involving  electrons  (hence  orbitals)  which  are  located  near  the 

molecule.  Because  of  this  constraint  the  models  involve  either  a 

19  20 

new  charge  transfer  (CT)  state  created  by  chemisorption,  '  or 
the  excitation  of  "metal  electrons"  into  empty  "molecular" 
orbitals . The  detection  of  a  CT  state  in  the  very 

37 

carefully  controlled  EELS  experiments  of  Avouris  and  Demuth 
gives  plausibility  to  the  idea  that  a  resonant  enhancement  of  the 
signal  can  play  a  role  in  SERS.  Unfortunately,  no  detailed  Raman 
measurement  has  been  made  on  that  system  at  a  frequency  that 
would  excite  the  CT  state.  It  is  conceivable  that  a  small 


transition  dipole  derivative  and/or  the  broadening  of  the  state 

12 

by  energy  transfer  to  the  surface  may  substantially  depress  the 

resonant  Raman  signal  and  result  in  a  small  enhancement. 

The  other  chemical  enhancement  mechanism  invoices  the  fact 

that  the  polarizability  of  the  metal-molecule  system  might  be 

much  larger  than  that  of  the  molecule  alone,  since  (qualitatively 

speaking)  some  of  the  metal's  electrons  which  can  scatter  from 

the  local  mode  have  much  higher  polarizability  than  the  electrons 

bound  to  the  molecule.  This  idea  is  at  the  basis  of  Otto's 
13  17 

proposal  '  for  the  role  of  atoms  in  SERS .  A  version  of  it, 
the  modulated  reflectance  model,  was  discussed  by  Maniv  and 
Metiu.^® 

In  the  present  paper  we  examine  theoretically  the  idea  that 
metal  atoms  on  a  solid  surface  might  increase  the  Raman 
polarizability  of  the  molecules  adsorbed  on  them.  To  estimate 
the  order  of  magnitude  of  this  increase  we  use  ab-initio  Hartree- 
Fock  calculations  for  a  CO  molecule  adsorbed  on  Cu.  This  method 
permits  us  to  do  a  better  job  in  computing  the. detailed  local 
binding  properties  than  the  Newns  -  Anderson  Hamiltonian  or  the 
jelllum  RPA  calculations  that  have  been  used  in  prior  work. 
Unfortunately,  the  large  number  of  electrons  severely  limits  our 
ability  to  deal  with  the  effects  of  the  extended  solid.  The 
following  strategy  was,  therefore,  employed.  Noting  that  atomic 
hydrogen  like  Cu  contains  a  single  S  electron  outside  of  a  core 
we  calculated  the  properties  (primarily  bonding,  but  also  the 

core-hole  states  of  H  CuCO  clusters  and  compared  them  with  the 
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results  of  Cu^CO  ab-initio  Hartree-Fock  calculations  ,  as  well 

as  with  the  experimental  results  for  CO  adsorbed  on  Cu  surfaces. 

It  turns  out  that  the  linear  H-Cu-CO  complex  reproduces  the 

properties  mentioned  above  reasonably  well.  Therefore  one  hopes 

that  the  polarizability  derivative  of  the  complex  with  respect  to 

the  C-0  stretch  will  also  agree  reasonably  well  with  experiments 

for  CO  bonded  to  a  Cu  atom  adsorbed  on  a  Cu  surface  (i.e.  the 

adatom  model).  We  do  not  expect,  of  course,  to  obtain  the 

correct  electromagnetic  properties  of  the  surface  bu^  this  is  not 


a  handicap  since  we  are  only  interested  in  the  local  modification 
of  the  polarizability  derivative  caused  by  the  chemical  bonding. 

Thus,  this  calculation  is  especially  relevant  to  Otto's  atom 

1  3,17 
model . 

zz 

The  main  result  is  that  Zee  where  z  is  along  the  CuCO 

axis  and  is  the  C-0  bond  distance,  is  increased  by  the 
binding  of  CO  to  Cu  by  a  factor  of  ^2.  This  gives  a  Raman 
intensity  enhancement  of  *-4 .  Careful  experiments^®  give  for  this 
quantity  an  order  of  magnitude  of  ~6,  or  of  less  than  10.  The 
experiments  were  not  carried  on  the  CO/Cu  system  and  we  mention 
them  here  only  to  indicate  that  we  are  in  general  agreement  with 
the  order  of  magnitude  obtained  from  data  taken  and  analyzed  with 
great  care . 

It  should  be  observed  that  without  adding  hydrogen  the  CuCO 
complex,  by  Itself,  gives  results  which  are  qualitatively 
different  from  those  of  CO  chemisorbed  in  an  on-top  site.  The 
hydrogren  atom  attracts  electrons  from  the  Cu  core  and  this 
allows  the  Cu  atom  to  bind  more  strongly  to  CO  as  it  does  in  the 
real  system.  By  carefully  choosing  the  cluster  it  may  be 
possible  to  accomplish  a  similar  simplification  for  studying 
chemisorption  properties  of  other  molecules  and/or  other  sites  on 
a  Cu  surface. 


II.  CLUSTER  MODEL 

31 

As  noted  in  the  Introduction  Bagus  and  Seel  (BS)  have 

previously  used  a  cluster  model  to  calculate  the  adsorption 

properties  of  CO  on  Cu(lOO).  They  considered  the  case  of  on-top 

binding,  which  is  appropriate  for  CO  in  a  C(2x2)  over layer 
32  33 

structure,  '  and  placed  the  Cu  atom  which  binds  to  CO  at  the 
apex  of  a  Cu^  square  pyramid.  The  four  Cu  atoms  in  the  base  of 
the  pyramid  represent  the  nearest  neighbors  in  the  second  layer 
of  the  (100)  copper  surface. 

Since  the  electronic  configuration  for  atomic  Cu  consists 
of  a  single  ^electron  outside  of  closed  shells  we  thought  it 
might  be  possible  to  simplify  the  BS  model  without  significant 
loss  of  accuracy  by  substituting  hydrogens  for  the  nearest 
neighbor  Cu  atoms.  However,  trial  calculations  on  CuH^  showed 
that  there  was  excessive  charge  transfer  (-0.56e)  from  the 
hydrogens  to  the  central  Cu.  Further  examination  revealed  (see 
Sec.  Ill)  that  this  charge  transfer  occurred,  in  large  part,  into 
the  4p  Cu  orbitals  of  tt  character  which  suggested  replacing  the 
four  atom  square  by  a  single  H  atom  on  the  symmetry  axis. 

Indeed,  upon  doing  so,  the  charge  transfer  was  reduced  to  '•O.OSe, 

a  reasonable  value  for  a  neutral  or  quasi-neutral  metal 

.  34 

surface. 

Of  course,  the  small  total  charge  transfer  is  just  one 

criterion  that  the  simplified  HCuCO  cluster  model  should  satisfy. 

It  also  ought  to  reproduce  other  properties  of  CO  binding  to  Cu, 

particularly  the  Cu-C  equilbrium  bond  length  and  the  Cu-C  bond 

energy.  As  we  will  see  in  Sec.  Ill  HCuCO  does  quite  well  for 

3  2  35 

both  of  these  properties.  In  fact,  compared  to  experiment  ' 
it  gives  a  better  binding  energy  than  the  BS  model  which 
underestimates  it. 

For  calculations  of  the  C-O  bond  polarizability  derivative 
it  is  important  that  the  valence  molecular  orbitals  primarily 
localized  on  CO  be  accurately  described.  Again,  we  will  see  in 
the  next  section  that  the  orbitals  obtained  from  our  model  are  in 
close  agreement  with  the  corresponding  Cu.CO  orbitals. 
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III.  COMPUTATIONS 

All  of  the  electronic  structure  calculations  reported  in 

this  paper  were  carried  out  by  the  ab  initio  Hartree-Fock  self 

36 

consisent  field  method  using  the  H0ND05  program.  For  open 

shell  systems  the  orbitals  were  unrestricted,  for. closed  shell 

systems  they  were  taken  to  be  doubly-occupied. 

A  suitable  CGTO  minimum  basis  set  for  describing  the  3^ 

transition  series  has  been  developed  by  Tatewaki  and 
37 

Huzinaga.  They  used  a  contraction  of  three  Gaussians  to 

represent  the  s_  and  g_  orbitals  but  found  that  a  fourth  GTO  was 

required  to  give  an  accurate  atomic  3^ orbital  energy.  This 

contracted  minimum  basis  will  be  referred  to  as  SET  A.  It  turns 

out,  however,  that  one  cannot  achieve  sufficient  binding  of  Cu  to 

CO  unless  a  4q_ polarization  function  is  added  to  the  basis.  To 

3  8 

this  end  we  adapted  the  4g,  orbital  obtained  by  Wachters  for  the 

3d^^4p  configuration  of  atomic  Cu.  Instead  of  eleven  Gaussian 

fvinctions  only  the  three  with  largest  coefficients  were  retained 

and  these  were  then  contracted  and  scaled  leading  to  basis  set  B. 

Xn  order  to  test  the  convergence  of  our  calculations  still  a 

third  set  (SET  C)  was  constructed  by  making  a  1/1/1  split  of  the 

4s_  orbital  and  a  3/1  split  of  the  3^. 

Three  different  contracted  Gaussian  bases  were  chosen  for 

the  CO  moiety  as  well.  The  first  {SET  1)  was  the  minimum  3G  set 
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of  Tavouktsoglou  and  Huzinaga  the  second  (SET  2)  was 

40 

Damning 's  double  zeta  [4s2p]  contraction,  and  the  third  a 

[5s4p]  set  (SET  3)  which  resulted  from  splitting  of  Dunning's 

41 

orbitals.  There  are,  all  together,  9  possible  combinations  of 
Cu  and  CO  basis  functions  which  may  be  uniquely  designated  A1 , 

...,  C3 .  Finally,  for  the  hydrogen  atom  we  employed  the  standard 

42 

Hehre,  Stewart  and  Pople  ST0-3G  basis  function. 

In  order  to  investigate  the  CO  basis  we  have  determined  a 
number  of  properties  for  the  free  CO  molecule  most  of  which  are 
shown  in  Tables  I  and  II.  From  the  tables  it  is  apparent  that, 
although  SET  1  is  inadequate,  SETS  2  and  3  reproduce  with  good 
accuracy  the  results  of  experiment  and/or  better  calculations  as 
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far  as  the  equilibrium  bond  length,  axial  polarizability 

Mulliken  gross  populations  are  concerned. 

The  error  in  the  dipole  moment  of  CO  compared  to  experiment 
43 

is  a  well-known  consequence  of  the  Hartree-Fock  approximation. 
Evidently,  the  consequences  are  not  as  serious  for  the 
polarizability  nor,  we  assume,  the  polarizability  derivative. 

44 

The  polarizability  was  obtained  by  finite  perturbation  theory 
using  a  quadratic  fit  to  the  energies  found  for  an  applied  field 
equal  to  -0.01,  0.00,  +0.01  a.u.  As  a  check  the  dipole  moment 
determined  by  fitting  agreed  with  the  directly  computed  value  (at 
F=0)  to  within  ±.002  a.u. 

By  comparing  basis  sets  2  and  3  in  Table  I  we  get  an 

indication  that  our  calculations  may  be  reasonably  well  converged 

with  respect  to  basis  set  size.  However,  it  is  hazardous  to  draw 

conclusions  in  this  regard  as  one  can  see  from  the  polarizability 
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resuts  of  McLean  and  Yoshimine^  and  Huo  quoted  in  the  table. 
Our  initial  test  of  the  Cu  basis  was  carried  out  on  the 

CuH.  cluster.  We  noticed  immediately  that  there  was  considerable 

4 

charge  transfer  from  the  hydrogens  to  Cu.  As  an  example,  the 

Mulliken  gross  populations  for  set  B  are  presented  in  Table  III. 

Clearly,  the  negative  charge  on  Cu  (0.56e)  is  much  too  large  for 

a  neutral  or  quasi-neutral  Cu  surface.  compared  to  the  d^^s^ 

free  atom  it  is  evident  that  the  additional  electronic  charge  has 

gone  into  the  4p  orbitals.  The  ir-like  4p '  s  ( 4p  ,  4p  )  alone  can 

X  y 

account  for  almost  the  entire  deviation  from  neutrality.  Thus, 
as  noted  earlier,  we  were  persuaded  to  substitute  a  single  H  atom 
on  the  symmetry  axis  for  the  H.  group.  Upon  doing  so  the  4p  and 
4Py  populations  fell  essentially  to  zero  which  was  expected  and 
the  charge  on  Cu  was  reduced  to  0.06  as  desired.  It  should  be 
observed  that  the  H  of  CuH  is  considerably  less  positive  than  the 

H's  of  CuH,  showing  that  the  effect  is  not  simply  due  to  having 

4 

only  1/4  as  many  H  atoms. 

One  might  wonder  whether  the  last  remaining  H  atom  can  be 
eliminated  as  well.  In  fact  it  must  be  retained  for  otherwise  CO 
does  not  bind  to  Cu,  at  least  in  calculations  with  the  basis  sets 


47 

used  here.  In  order  to  see  whether  the  HCuCO  model  adequately 

remedies  this  deficiency  we  did  preliminary  computations 

employing  the  A1  and  B1  bases.  Only  the  Cu-C  distance  was 

varied,  the  C-0  and  Cu-H  bond  lengths  were  fixed  at  the 

48 

equilibrium  minimum  basis  set  value  for  the  free  molecules. 

The  binding  parameters  that  were  found  are  compared  in  Table  IV 
with  those  of  BS  and  experiment.  Without  4p  orbitals  in  the  Cu 
basis  (SET  A1 )  the  binding  is  too  weak.  Once  such  polarization 
functions  are  included  (SET  Bl)  the  results  are  quite 
satisfactory  although  the  excellent  agreement  with  experiment  may 
be  partly  fortuitous.  A  Mulliken  gross  population  analysis  of 
the  three  highest  valence  orbitals  derived  from  free  CO  (Of. 

Table  V)  -  carried  out  for  SET  B2  -  further  demonstrates  the 
similarity  between  the  BS  model  and  ours  as  far  as  the  CuCO 
entity  is  concerned. 


IV.  RESULTS  AND  DISCUSSION 

Our  final  results  for  the  HCuCO  cluster  model  of  CO 

adsorbed  on  Cu(lOO)  are  reported  in  Tables  VI  and  VII.  The 

former  table  pertains  particularly  to  the  properties  of  the  Cu-C 

bond,  the  latter  to  the  C-0  bond.  In  both  cases  the  B2  and  C2 

basis  sets  were  employed  along  with  a  Cu-H  bond  length  fixed  at 

2.709  a.u.  The  latter  distance  is  the  equilibrium  value  for  the 

free  molecule  in  the  minimum  basis.  This  is  close  to  the 

49 

experimental  result  of  2.765  a.u.  and  we  would  not  expect  it  to 
change  much  with  enlarging  the  basis  set  or  forming  the  HCuCO 
cluster.  We  also  tested  this  calculation  to  insure  that  the 
polarizability  derivative  we  seek,  6a^^/6R_  is  insensitive  to 
small  variations  in 

In  Compiling  Table  VI  the  C-0  bond  length  was  held  constant 
at  the  free  molecule  value  determined  by  the  SET  2  basis.  Note 
that  the  binding  properties,  which  are  similar  in  value  to  our 
preliminary  results  (Cf.  Table  IV)  remain  in  good  agreement  with 
experiment.  The  Cu-C  bond  polarizability  derivative  is  quite 
small  as  may  be  seen  by  comparing  the  derivative  for  the  C-0 
bond.  We  used  the  Cu-C  equilibrium  bond  distances  computed  for 
this  table  as  input  to  obtain  the  C-0  bond  properties  given  in 
-  Table  VII. 

a_  is  a  nearly  linear  function  of  R_„  _  and  R_  _  in  the 
vicinity  of  the  equilibrium  bond  lengths.  Hence,  the  specific 
values  chosen  for  these  parameters  is  not  important.  If  is 

increased  from  2.136  to  2.256  a.u.,  for  example,  the  (SET  C2) 
polarizability  derivative  in  Table  VII  changes  from  16.51  to 
16.80  a.u. 

The  differences  between  basis  sets  B2  and  C2  again  provide 
some  indication  of  how  well -converged  the  results  are.  In 
general,  these  differences  are  relatively  small.  Although  the 
latter  statement  is  not  true  for  the  polarizability  derivative 
6a  /6R  in  that  case  the  change  is  small  on  an  absolute 

ZZ 

basis. 

'As  we  have  already  emphasized  in  the  Introduction/ our' 


intention  was  to  estimate  to  what  extent  the  non-resonant  Raman 
signal  for  the  C-0  stretch  is  increased  when  CO  binds  to  a  Cu 
surface.  An  HCUCO  cluster  has  been  used  to  simulate  the  on-top 
binding  of  CO  on  Cu(lOO).  The  fact  that  our  cluster  reproduces 
the  binding  properties  of  adsorbed  CO  and  the  CO  valence  orbitals 
of  pyramidal  Cu^CO  gives  us  hope  that  a  correct  estimate  of  the 
Raman  intensity  has  been  obtained.  This  is  especially  likely 
since  the  property  sought  is  a  local  quantity.  Our  result  is 
that  the  'Chemical '^enhancement  is  about  a  factor  of  four  which 
is  in  general  agreement  with  the  few  careful  estimates  based  on 
experimental  data.  '  We  emphasize  that  the  present 
calculation  does  not  deal  with  the  question  of  resonant 
enhancement.^ 
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Mulliken  gross  populations  for  the  three  highest 
valence  molecular  orbitals  of  free  CO  as  a  function 
of  basis  set. 


CARBON  OXYGEN 


Table  III.  Mullilcen  gross  populations  in  CuH  and  CuH  obtained 
with  the  SET  B  basis. 
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FIGURE  CAPTIONS 

Schematic  diagram  for  Raman  scattering  by  a  diatomic 
adsorbed  on  a  metal  surface.  The  shape  of  the 
electron  cloud  is  drawn  approximately.  Lined 
circles  represent  ion  cores;  wavy  lines  photon 
propagators;  full  lines  are  electron  propagators, 
and  the  double  line  is  a  phonon  (i.e.,  the 
stretching  mode  of  the  adsorbed  diatomic  molecule 
propagator) .  The  energies  for  each  propagator  are 
noted  along  the  lines.  A  more  detailed  description 
of  the  process  is  given  in  the  text. 


